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Abstract 

The Weibel instability could be responsible for the generation of magnetic fields in various objects 



o 
o 

such as gamma-ray bursts, jets from active galactic nuclei, and clusters of galaxies. Using numerical 

O 

simulations, the development of the Weibel instability at a temperature gradient is studied. It is 

CO 

found that current sheets are first generated at the gradient, and then they are rounded off and 



turn into current filaments. During this process, return currents are generated around the filaments 



and they prevent filaments from merger. The magnetic fields around the filaments persist at least 
until t ~ 8000/ 'u p , where u p is the plasma frequency, and it is very likely that they survive for a 
much longer time. 
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I. INTRODUCTION 



The question of the origin of magnetic fields in the Universe is one of the most challenging 
problems in modern astrophysics. One of the fascinating ideas is that the fields are generated 
by the Weibel instability, a plasma instability in a collisionless plasma [lj. This instability 
is driven by the anisotropy of the particle velocity distribution function (PDF) of plasma. 
When the PDF is anisotropic, currents and then magnetic fields are generated in the plasma 
so that the plasma particles are deflected and the PDF becomes isotropic |2[. Through 
this process, the free energy attributed to the PDF anisotropy is transferred to magnetic 
field energy. This instability does not need seed magnetic fields. It can be saturated only 
by nonlinear effects, and thus the magnetic fields can be amplified to very high values. 
In particular, the generation of magnetic fields at shocks through the instability has been 
studied by many authors, because the PDF is anisotropic at shocks 

However, the magnetic fields generated by the Weibel instability at a shock may not 
survive for a long time. At a shock, current filaments are created and magnetic fields are 
generated around the m flQ. As the — mff e together by the magn et 1C force between 
them, the magnetic field strength increases j?l llj However, the magnetic field should be 
saturated when the current strength reaches the Alfven current llll . which is the maximum 

n 

current allowed by the self-generated magnetic field [12j|. The Alfven current is given by 

I A = (mc 3 /q)j(3 , (1) 
where m and q are the mass and charge of a particle, respectively, f3 is the mean velocity 



of particles normalized by the speed of light c, and 7 = 1/y 1 — [3 2 . After the saturation, 
magnetic field strength, B, would decrease as the filament size r increases through mergers, 
because B oc Ia/t and I a is constant, although current numerical simulations cannot fully 
deal with this phase because of the limitation of the simulation box size. The timescale of 
he filament mergers before the saturation is ~ 10 w" 1 , where uo p is the plasma frequency 
lol |. If the timescale after the saturation is the same as that before the saturation, the 
magnetic fields would rapidly fade away. 

Another site where the Weibel instability could be effective is temperature gradients, 
where the PDF is also anisotropic jl^. One example of such temperature gradients is cold 
fronts observed in clusters of galaxies [3, [ljj]. In this paper, we present the results of 
numerical simulations of electron-positron plasma performed to investigate the long-term 



evolution of the Weibel instability at a temperature gradient. We emphasize that in our 
simulations, Coulomb collisions are ineffective (coUisionless plasma), which is different from 
the assumption of Ref . [k| . 

II. MODELS 

We performed numerical simulations of an electron-positron plasma at a temperature 
gradient. The simulation code used is a relativistic, electromagnetic, particle- in-cell code 
with two spatial and three velocity dimensions, which was developed based on a general 
description by Ref. 151 ] . The code is a momentum conserving code. Using the code, we 



solve the Maxwell equations (in Gaussian units): 



1 dE Aty 



V • E = Anp , V • B = , (4) 

where E is the magnetic field, B is the magnetic field, J is the current density, and p is the 
charge density. We also solve the equation of motion for each particle: 

* J E + »±*), (5) 



dt \ •ym e c 

where p is the momentum of a particle, and m e is the electron mass. The simulations are 
performed on a 1024 x 512 grid (the axes are labeled as x and y, respectively) with a total 
of 20 million particles. Temporal and spatial scales in the simulations are normalized to 
the inverse electron plasma frequency uj~ 1 = (47rn e oe 2 /m e )~ 1//2 and the coUisionless skin 
depth A e = c/ujp, where n e0 is the average initial electron or positron density, — e is the 
electron charge, and c is the speed of light, which is the unit of velocity in our simulation 
code. The units of mass and charge are the electron mass m e and the absolute value of 
the electron charge e, respectively. In these normalized units, the box size is 160 x 80. For 
electromagnetic fields, we adopt a periodic boundary condition. For the x direction, we set 
walls at x = 20 and 140. The regions, < x < 20 and 140 < x < 160 are used for Joule 
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dissipation of electromagnetic waves so that electromagnetic waves do not enter from the 
other side. 

For particles, we adopt a periodic boundary condition for the y direction. On the other 
hand, for the x direction, we adopt reflection boundary conditions with 'heat walls'. Particles 
that hit the wall at x = 20 have a thermal velocity of = 0.1. That is, the reflected particles 
have the Maxwellian velocity distribution with the deviation of 0.1. At the wall of x ~ 140, 
reflected particles have the one with the deviation of ctr = 0.5. For the right wall, we 
randomly choose the reflection point in 137 < x < 140 to avoid the formation of artificial 
structures in the plasma at the wall. Because of these boundary conditions, the total energy 
of the system does not conserve. 

Initially, the plasma in the simulation box has no magnetic fields and is 'isobaric', that is, 
the square of the thermal velocity ('temperature') times the density is constant. At the left 
(x = 20) and right (x = 140) boundaries, the plasma has thermal velocities of = 0.1 and 
or = 0.5, respectively. For 20 < x < 140, the temperature changes linearly as x increases. 
Thus, the density on the left side is higher. 

It is to be noted that in another simulation we have calculated the plasma evolution when 
two plasmas with different temperatures but with the same pressure are bordered at x — 80. 
The results are qualitatively the same as those shown below. 

III. RESULTS 

Fig. ^ shows the evolution of the magnetic energy, Wb, in units of the total energy in 
the box at t = (the sum of the kinetic energies of all particles). For comparison, the total 
particle energy, Ep, is shown. The total particle energy decreases because of the cooling at 
the left wall as well as the generation of magnetic fields. 

The magnetic energy rapidly increases at t < 400. Because of the initial density gradient, 
particles as a whole move in the right direction. In particlar, those with large velocities 
move fast. Thus, on the frame moving with these particles, the effective temperature of the 
particles is smaller in the x-direction than those of the other directions, which develops the 
Weibel instability. This initial stage of the instability will be studied in detail elsewhere. 
Fig. 121 shows the current density in the ^-direction (J z ) at t = 400. Current sheets are seen 
at x < 80. Strong magnetic fields are formed around these current sheets. For x > 80, 
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FIG. 1: The evolution of the magnetic energy (Wb' solid) and the total particle energy (Ep: 
dotted) in units of the total energy in the box at t = 0. 
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FIG. 2: (Color online). The current density in the z-direction ( J z ) at t = 400 in units of n e Q ec. 

sheets are still developing. On the other hand, the current density in the x-direction (J x ) is 
much weaker than J z and does not contribute to the formation of magnetic fields. 

Note that contrary to the Weibel instability at the temperature gradient, current sheets 
are not created at shocks. While the thermal velocity of particles is larger in the 'directions' 
perpendicular to the temperature gradient in the calculation presented here, it is larger in 
the 'direction' of the shock normal at a shock front Q|. This is the reason why only current 
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FIG. 3: The same as Fig. |2 but for t = 5000. 

sheets (two-dimensional structures) are generated in the former. In the latter, only current 
filaments are formed in the direction of the shock normal. 

At t ~ 400, Wb reaches its local maximum. We found that the thickness of individual 
current sheets at x < 60 in Fig. El is comparable to the gyroradius of particles. Note that for 
current filaments, if the radius of the individual filaments is comparable to the gyroradius 
of particles, the current strength is to be the Alfven current ^l]. It has been shown that 
the magnetic field strength should reach its maximum at that time 

At 650 < t < 1100, the mixing of hot and cold plasmas proceeds, and Wb decreases. 
However, for t > 1100, Wb starts increasing again. In this period, the current sheets are 
rounding and turning into 'filaments', because the sheets are unstable. Since the sheets are 
perpendicular to the temperature gradient, the filaments should also be perpendicular to 
that direction. In this process, the cross sections of the sheets decrease. The magnetic field 
strength around the filaments slowly increases through the rounding and shrinking. The 
increase of the magnetic fields produces electromotive forces around the filaments, which 
induce return currents around the filaments. The return currents are clearly seen around 
filaments at t > 1500 (Fig.©. 

Since the return currents shield the magnetic fields outside the filaments, the magnetic 
force between the filaments is significantly reduced. At this stage, the filaments stop merging. 
While the filaments at large x do not evolve, the sheets continue to round at smaller x, and 
Wb gradually increases even at t > 2000 (Fig. [T]). Fig. 01 also shows that current density is 
not uniform in a filament. The current density is small at the center of the filament and 
thus each filament appears to be 'a current tube'. 

In Fig. El filaments are seen only at x > 80. The thermal velocity outside the filaments 
does not much depend on x and is ~ ol — 0.1. Particles with large velocities come from the 
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FIG. 4: (Color online). The number density of electrons at t = 5000 in units of n e o- 

right wall. However, the magnetic field around the filaments prevents those particles from 
moving to the left wall. Thus, the PDF anisotropy is larger at larger x, which is the reason 
that filaments develop at larger x. 

In our simulations, the particles are merely reflected and redistributed thermally at the 
right and left walls; they are not uniformally fed into the system. In spite of this, the walls do 
not create artifacts. This may be because particles in low and high density regions bordered 
on the same wall are well-separated by magnetic fields and they are not mixed together. 
It would be interesting to compare the results here with those of the simulations in which 
particles are uniformally fed into the system. 

IV. DISCUSSION 

We studied the structure of the stable filaments seen at x > 100 in Fig. El In the 
filaments close to the right wall, the average drift velocity of particles in the ^-direction is 
|An,z| ~ 0.05, while their thermal velocity is o m ~ 0.16. This means that the filaments 
are 'hot beams'. Figs. 0] and show the number density and the kinetic energy density of 
particles at t — 5000, respectively. For the filaments at x > 100, while the number density is 
small inside them, the kinetic energy density is almost the same between the inside and the 
outside of the filaments. Fig. |H1 shows the magnetic energy density (Ub), and demonstrates 
that magnetic fields are generated only at the surfaces of the filaments. There is a relation of 
Ub = (1/2)(B / B*) 2 , where = cy/47m e om e . The magnetic field strength is determined so 
that the gyroradius of particles with the velocity of a m ~ 0.16 is comparable to the thickness 
of the current tubes. This magnetic field confines hot particles within the tubes. 

Based on these facts, we can analyze the structure of the magnetic fields as follows. We 
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FIG. 5: (Color online). The kinetic energy density of particles (electrons and positrons) at t = 5000 
in units of n e om e c 2 . 
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FIG. 6: (Color online). The magnetic energy density, Ub, at t = 5000 in units of n e om e c 2 



use the Gaussian units here. We define 'thermal pressure' as 



P = 2n P m P a 2 c 2 



(6) 



where n e is the density, and a is the thermal velocity of particles. The factor of two comes 
from the fact that there are electrons and positrons. Our simulations showed that the 
thermal velocities inside and outside a filament are a in < a R and cr out > (a in > cr out ), 
respectively. Fig. suggests that 

P m m P out , (7) 

where P in and P out are the thermal pressures inside and outside a filament, respectively. 
Equation (jZJ) means that 



n 'm°']n ~ n ont cr , 



out J 



(8) 



where n in and n out are the electron densities inside and outside a filament, respectively. 
Since the gyroradius of particles inside a filament is comparable to the thickness of the 
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current tube (A), the magnetic field strength is given by 

2 

m e c <T[ n . . 

B « V— z , (9) 

Ae 

where rj is the ratio of the thickness to the gyroradius. On the other hand, the magnetic 
field strength is also given by B = 21 /(re), where I is the current in the filament, and r 
is the radius of the filament (or the tube). If we assume that the average drift velocity of 
particles in a filament is given by |An,z| = £cin, the current is I — 47m in e/3i niZ crA. Thus, the 
magnetic field strength is 

B = 8n£n in ea in A . (10) 
Combining equations (JSJ), (JHJ), an d (fTUj) . we obtain 

B « y/tojiB^t , (11) 

and 

A^ A /^A e ^. (12) 

V 2^ V out 1 ; 

We used the relation n out m n e0 . In our simulations, we observed <jj n — 0.16, cr ut — 

0.08, and 

£ = 0.3 for the filaments close to the right wall. For these values and r) — 2, equations (fTTj) 
and (|T2*|) predict that B/B* m 0.09 and A/A e ~ 4, which are consistent with the results 
of the simulations (Fig. The above equations alone cannot determine the radius of a 
filament r. It may depend on the initial conditions of plasma. 

Although we finished the simulations at t ~ 8000, we expect that Wb continues to increase 
for t > 8000, because of the rounding and shrinking of sheets at smaller x. On the other 
hand, the filaments at larger x will remain the same. Fig. shows the histogram of 7 — 1 
(the kinetic energy of particles in the unit of m e c 2 ) at t = 8000. It shows that particles 
consist of those with the thermal velocity of ~ &l and those of ~ gr. A sign of particle 
acceleration cannot be found. 

Our simulations are two-dimensional in space. Thus, we cannot deal with three- 
dimensional deformations of current filaments, which might be caused by the Kink instability 
or filament mergers. However, for the mergers, if the time-scale of the formation of return 
currents is smaller than that of the mergers, the results presented here would not much 
change. For filaments that have not been shielded by return currents, it would take a longer 
time to merge together, because the filaments would not be aligned and the magnetic inter- 
action among them would be weaker in a three-dimensional space. Thus, filament mergers 
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FIG. 7: The histogram of 7 — 1 (the kinetic energy of particles in the unit of m e c 2 ) at t = 8000. 



would be further prohibited. On the other hand, while our simulations are two-dimensional 
in space, they are three-dimensional in velocity space. Therefore, the resultant velocity 
distribution of particles would not be much different from that in full three-dimensional 
simulations. 

We have shown that magnetic fields generated at a temperature gradient could survive 
as long as the gradient exists. The essence is the generation of current sheets, which later 
round and turn into filaments. Return currents created in this process shield the filaments 
and prevent them from mergers. It would be interesting to study this instability for electron- 
proton plasma. Assuming that protons and electrons have the same temperature, thermal 
velocity of the protons is smaller than that of the electrons by a factor of ^m e /m p , where 
■trip is the proton mass. Since for protons is larger than that for electrons by a factor of 
^Jm p /m ei the resultant magnetic fields would not be much different from that for electron- 
positron plasma (see Eq. [11]). although it must be confirmed by numerical simulations. 
It would also be interesting to apply the results to actual objects in the Universe, such as 
gamma-ray bursts, supernova remnants, and clusters of galaxies. However, the temperature 
gradients in them are much smaller than that studied in this paper, which might make some 
difference. Moreover, the coherent scale of the magnetic fields generated by the Weibel 
instability is much smaller than astrophysical scales. Thus, some mechanisms that increase 
the scale would be required Q. In the future, we would like to challenge these problems. 
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